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Nearly degenerate four-wave mixing (NDFWM) within a 
closed degenerate two-level atomic transition is theoretically 
and experimentally examined. Using the model presented by 
A. Lezama et al [Phys. Rev. A 61, 013801 (2000)] the ND- 
FWM spectra corresponding to different pump and probe po- 
larization cases are calculated and discussed. The calculated 
spectra are compared to the observation of NDFWM within 
the 6S'i/2(-F = 4) ^ (iP3/2{F = 5) transition of cesium in a 
phase conjugation experiment using magneto optically cooled 
atoms. 



I. INTRODUCTION. 

The generation of new optical fields is the most spec- 
tacular consequence of the intrinsic non-linearity of the 
light-matter interaction. Four- wave mixing (FWM), i.e. 
the generation of a fourth field as the result of the inter- 
action of a material sample with three electromagnetic 
fields is the simplest non-linear optics generation process 
allowed in every material . 

This paper is concerned with nearly degenerate four 
wave mixing (NDFWM) occurring when two of the three 
incident fields (designated as pump fields) have the same 
frequency loi and the third field (probe field) has an inde- 
pendent frequency u}2 = uji + 5 which is tunable around 
bJi. Of particular interest is the case where all fields are 
nearly resonant with the same optical atomic transition. 
In such case the spectroscopic study of the NDFWM sig- 
nal as a function of the probe to pump frequency offset 
can provide useful information about the atomic dynam- 
ics. NDFWM was used by Rothberg and Bloembergen ||] 
for studying coUisional dynamics in sodium vapor. Boyd 
et al Q] have studied theoretically the atomic polarization 
and the propagation of the probe and NDFWM waves in 
the case of a pure two-level system (PTLS) driven by a 
arbitrarily intense pump. Their work illustrate the crit- 
ical influence of the atomic damping mechanism in the 
FWM process. Steel et al [|| explored NDFWM in both 
open and closed PTLS. Berman and coworkers Q clearly 
established the suitability of NDFWM for the determi- 
nation of all relevant relaxation parameters. Andersen 
and coworkers analyzed the NDFWM process in PTLS 
in terms of the dressed state picture 0. Their work 



stresses the importance of the initial atom-field state on 
the FWM yield. 

In spite of the fact that the experiments on NDFWM 
generally involve degenerate energy levels, most experi- 
ments to date were carried under conditions where the 
level degeneracy does not play an essential role. In 
such cases, simplified theoretical models using two-level 
schemes were successfully used. Nevertheless, the role of 
the level degeneracy for NDFWM was clearly appreci- 
ated by Berman and coworkers Q who stressed the im- 
portance of the non-conservation in the evolution of the 
multipolar expansion terms of the atomic density ma- 
trix, for the observation of sub-natural resonances in the 
NDFWM signal spectrum. Their work include a detailed 
calculation, valid for degenerate atomic levels, of the non- 
linear atomic polarization to the lowest order in the in- 
cident fields. Recently, a model allowing the numerical 
investigation of the complete response of a degenerate 
two-level system (DTLS) to first order in the probe field 
for an arbitrarily intense pump and arbitrary pump and 
probe polarizations was presented. 

Recent experimental observations illustrate the need 
for a deeper understanding of the role of the level de- 
generacy in the non linear response of DTLS. An impor- 
tant example is provided by the spectroscopic response 
of alkaline atoms in a magneto-optical trap (MOT). In 
this system the atoms are submitted to the excitation 
of the trapping beams which are quasi resonant with a 
closed degenerate transition. The absorption [^,0 and 
FWM |]lT| , |l2| spectra of magneto-optically trapped atoms 
present characteristic features resulting from the level de- 
generacy. Recently, the use of DTLS for efficient and 
highly selective NDFWM generation in an atomic vapor 
was reported jl^ . 

The aim of this paper is to address the role of the level 
degeneracy in NDFWM on a closed transition between 
two degenerate atomic levels. The dependence of ND- 
FWM generation on the intensity of the pumping field 
and the polarizations of the pump and probe will be an- 
alyzed. In the next section, NDFWM spectra calculated 
for an homogeneous ensemble of DTLS under different 
conditions are presented. The following section is de- 
voted to the experimental observation of NDFWM in cold 
cesium atoms. Discussion of the experimental results in 
regard to the theoretical predictions follows. 
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II. THEORETICAL PREDICTIONS. 

The NDFWM generation in an homogeneous ensemble 
of DTLS is analyzed using the semiclassical procedure al- 
ready presented in In Fig. |l| the essential elements 
of this model are reminded. The two-level atoms are 
driven by a pump field of frequency uji of arbitrary in- 
tensity and polarization. The reduced pump field Rabi 
frequency is fli. The driven atoms are probed by a weak 
field of frequency lu2 = lui + 5 and arbitrary polariza- 
tion. Both fields are nearly resonant with the atomic 
transition (transition frequency IIuja) between the ground 
level g of total angular momentum Fg and the excited 
level e of total angular momentum fg. Having in mind 
the experimental observation presented below, we will 
restrict our attention to closed atomic transitions with 
Fe = Fg + 1. Level e decays back into g through sponta- 
neous emission at rate F. In addition to this relaxation 
mechanism a state independent decay rate 7 (7 ^ F) 
is assumed which represents the departure of the atoms 
from the interaction zone. This departure is compen- 
sated at steady state by the arrival (pumping term jpo) 
of "fresh" atoms isotropically distributed in the ground 
level, (po is the system density matrix in the absence of 
any applied field). Since no other relaxation mechanism 
exist for the ground state, 7 effectively plays the role of 
a ground state relaxation rate. The calculation is exact 
to all orders in the pump field and to first order in the 
probe field. The spectra represent the variation of the 
squared modulus of the atomic polarization induced at 
the frequency lo^ = 1lo\ —lu2 = 1^1 — ^ as a function of the 
probe to pump frequency offset S. Since the calculation 
corresponds to the response of an homogeneous ensem- 
ble of atoms at rest, no propagation effects such as phase 
matching, spatial emission pattern of the generated wave 
or beam absorption (or amplification) are accounted for. 

Before presenting the spectra calculated for DTLS, let 
us remind the main features of the NDFWM spectra for 
PTLS Figure | represent the spectra of the ND- 

FWM power as a function of ^ for A = and A = 2F 
(A = u!A — wi) for different values of the pump field 
Rabi frequency fio. We restrict ourselves to closed tran- 
sition where spontaneous emission is the only relaxation 
mechanism. For all values of f2oj the spectra present two 
symmetric sidebands at S = ztft where il = y^fl^ + A"^ 
is the generalized Rabi frequency. The width of these 
two peaks is determined by F. When fio ^ A,F a third 
peak develops at 5 = whose width is also determined by 
F. Only when O ~ F the central peak overcomes the two 
sidebands. As discussed in , the relative amplitude of 
peak at 5 = is determined by the coherence relaxation 
rate. In the present case this rate equals F/2 since only 
radiative decay is assumed. When the transition is open 
and when the escaping rates out of the two-level systems 
are different for the upper and lower level, an additional 



peak appears at ^ = with a width determined by the 
lower level relaxation rate . Finally let us mention that 
in PTLS the maximum NDFWM power is obtained for 

We examine now the NDFWM power spectra for a 
DTLS with Fg ^ 1 and F^ ~ 2 (no magnetic field 
present). Three polarization cases were considered: cir- 
cular and equal pump and probe polarizations, paral- 
lel linear pump and probe polarizations and linear and 
perpendicular pump and probe polarizations. The cor- 
responding level schemes are shown in Fig. |^ where an 
appropriate choice of the quantization axis is made in 
each case. 

Figure ^ shows the spectra calculated for pump and 
probe fields with the same circular polarization (cr^). As 
expected, due to the optical pumping of the population 
towards the Zeeman sublevel in g and e with the highest 
magnetic quantum number (see Fig. ^ a), the system 
behaves as a PTLS in the limit of large Oi. However, 
striking differences appear for fii < F in which case the 
Zeeman optical pumping is only partial. Notice the pres- 
ence at (5 = of narrow resonances: a peak for fti = F, 
A = 2F and a dip for 17i = F, A = and f^i = 2r, 
A = 2F. The width of this narrow resonance is deter- 
mined by 7. In this configuration, the pump and probe 
fields interact with pairs of Zeeman sublevels of e and g 
with me = Trig + 1 (rrii is a magnetic quantum number). 
Coupling between different pairs occurs only via sponta- 
neous emission. NDFWM signal arises from the coherent 
contribution of different Zeeman sublevels pairs excited 
by the field. Due to spontaneous emission, a given pair 
is not a close system in which case, a narrow structure 
(width 7) is expected at ^ = [^. As observed in Fig. 
I the sign of this narrow feature may vary depending on 
the pump field intensity and detuning. It is the result of 
the quantum interference of the contributions from the 
different Zeeman sublevels pairs to the non linear atomic 
polarization. A similar situation occurs in the case of 
parallel linear polarizations discussed next. 

Figure |^ refers to the spectra calculated for pump and 
probe fields with parallel linear polarizations. Unlike the 
previous case, in the limit of large fii, the system does 
not behave as a PTLS. Instead, taking as the quantiza- 
tion axis the direction of the common pump and probe 
polarization, it should be seen as three (2_Fg + l) two-level 
systems linked through spontaneous emission (Fig. ^ b). 
Each two-level system consist of a two Zeeman sublevels 
with rUg = rrie. Since the strength of the atom field cou- 
pling dependents only on |mg|, two different values of the 
pump Rabi frequency occur in the present example. This 
explains the splitting observed in the two sidebands for 
large fli. At low pump intensities a narrow resonance 
(width determined by 7) appears at (5 = for A = 2F. 
Also a narrow dip is present at 6 — for A = and 
fii < F (not shown in the figure). The appearance of 
narrow features for zero pump to probe detuning with 
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a width determined by the effective ground state relax- 
ation rate was predicted for open two- level systems [^. 
In the present case this feature should be seen as the 
consequence of the population transfer, due to sponta- 
neous emission, between pairs of levels connected by the 
fields. A detailed look into the evolution of this narrow 
resonance with the pump field power (for A = 2r) is 
presented in Fig. ^. The resonance at (5 = is nega- 
tive (dip) for fli > 0.4r. This result is to be compared 
with the contribution to the nonlinear atomic polariza- 
tion arising from the different pairs of Zeeman sublevels 
with rrie — rrig. These contributions can be extracted 
from the calculation after identification of the density 
matrix coefficients corresponding to a given value of \mi\ 
[i = e,g). Figs. and ^ represent the square modulus 
of the nonlinear atomic polarization contributions aris- 



ing from the Zeeman sublevels pairs with 



and 



\mg\ = 1 respectively. Only in the case of \mg\ = and 
for the largest value of Sli the central resonance is op- 
posite to the lateral sidebands. In all other cases this 
resonance has the same sign than the lateral sidebands. 
This is always the case for \mg\ = 1. Consequently, the 
dip observed at (5 = in the NDFWM spectra in Fig. ^ is 
the result of (destructive) quantum interference between 
the two contributions coming from \mg\ =0,1. Indeed, 
the two contributions have opposite phase over all the 
considered range of ili. 

The NDFWM spectra obtained for perpendicular lin- 
ear pump and probe polarizations (Fig. ^) are rather 
different to the ones presented above . For all the con- 
sidered pump field intensity range the spectra are domi- 
nated by the features occurring around 6 = 0. In fact for 
ill ^ 2r the spectral sidebands are barely visible and the 
spectra is mainly composed by a large central resonance 
with a width of the order of 7 |^ . This configuration pro- 
duces for ni^T the largest NDFWM yield For in- 
creasing fli the central peak splits into two components. 
Also in this case the spectra are asymmetric for A 7^ 0. 
The main features of these spectra can be qualitatively 
understood using the dressed states picture of the degen- 
erate atomic system in the presence of the pump field. A 
similar analysis was carried in [p^ (Appendix) to examine 
the probe absorption spectra of driven DTLS in this con- 
figuration. The dressed-state energy level scheme for a 
Fg = 1 ^ Fe = 2 transition driven by a tt polarized pump 
is presented in Fig. |l^ following the conventions adopted 
in jl^ . NDFWM resonances are expected to occur when 
the probe field is resonant with a transitions between 
dressed levels that are coupled to the pump field pho- 
tons. Four different probe frequencies satisfy this condi- 



tion. The corresponding values of S (see Table II in |15 



are indicated with solid arrows in the two lower spectra 
of Fig. ^ They correspond to the main features in the 
NDFWM spectra. From a simple dressed-state analysis, 
no resonance is expected to occur for the probe frequen- 
cies corresponding to transitions ending in the 1 1" , TV) 



dressed-atom levels since these states are not coupled to 
the pump photons due to the nie = nig selection rule 
valid for a tt polarized pump. Nevertheless small fea- 
tures appear in the calculated spectra at these positions 
(dashed arrows in Fig. ^). They are due to non-secular 
terms usually neglected in the dressed atom approach. 



III. EXPERIMENTAL OBSERVATION OF 
NDFWM IN COLD CESIUM. 

The observation of the NDFWM spectra was per- 
formed in a sample of cold cesium atoms produced in 
a magneto-optical trap (MOT). Light from a Ti:sapphire 
laser, nearly resonant with the cesium cycling transition 
65'i/2(i^ = 4) ^ 6P3/2{F = 5), was employed both for 
the trapping of the atoms and for the investigation of 
the NDFWM. The experimental setup is shown schemat- 
ically in Fig. The frequency of the Ti:sapphire laser 
is red-detuned by approximately two natural linewidth 
(r/27r = 5.3 MHz). A repumping diode laser, not 
shown in the figure, recycles the population lost to the 
65'i/2(F = 3) ground state. We use a backwards FWM 
configuration, where the two counter-propagating pump 
fields, the forward (F) and backward (B), have the same 
frequency and the same linear polarization. The probe 
beam (P), is linearly polarized and makes a small angle 
(0 — 4°) with the pumping beams. The probe beam has 
its frequency scanned around the frequency of the pump 
beams with the help of two acousto-optic modulators as 
shown in Fig. The relative polarization between the 
pump and the probe fields is controlled by a half-wave 
plate. The generated (nearly) phase conjugated signal 
(PC), which propagates in opposite direction with re- 
spect to beam P, is refiected out of a 50/50 beam splitter 
and detected by a fast photodiode. The trapping beams 
are switched on and off by a mechanical chopper with 
a transmission duty cycle of 95%. The number of cold 
atoms was estimated by measuring the absorption of the 
probe beam and is of order of lO"^. The NDFWM spectra 
were recorded within a ~ 1 ms time interval during which 
the trapping beams were blocked and the quadrupole 
magnetic field was turned off. Each of the F and B pump 
beams have an intensity of 7 mW/ cm? . The probe beam 
intensity is approximately equal to 0.7 raWlcm? . No sig- 
nificant modification of the spectra was observed at lower 
probe beam intensities. The maximum NDFWM power 
generated was of the order of 1 fiW . Typical spectra, 
recorded as a function of the pump to probe frequency 
offset (S, are shown in Figs, [l^a) and [l^a) for the probe 
beam polarization respectively parallel and perpendicu- 
lar to the pump beam polarization. 
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IV. DISCUSSION. 

The comparison of the observed spectra with the the- 
oretical predictions presented above is not direct. While 
the calculation applies to a homogeneous ensemble of 
atoms at rest, the PC signal is the result of coopera- 
tive emission of atoms under different excitation condi- 
tions. In the experiment the pump field, at a given po- 
sition in the sample, is due to the combined incidence of 
the F and B beams. The counterpropagating geometry 
of these beams produces a standing wave, consequently, 
the pump field intensity is spatially modulated. In gen- 
eral the pump field polarization may be spatial depen- 
dent. However, in the present study we have restricted 
ourselves to the case where the beams F and B have 
the same linear polarization and consequently no spatial 
variation of the polarization occurs. Also, the PC signal 
is sensitive to propagation effects in the atomic sample 
which are responsible for the phase matching condition 
and for possible spatial dependent amplification or de- 
pletion of the NDFWM through the medium. Finally, 
one should remind that the phase matching condition 
imposes constraints on the polarization components that 
may be present in the PC beam (only transverse compo- 
nents). However, in the two polarization cases considered 
in the experiment, linear parallel and linear perpendic- 
ular pump and probe polarizations, due to symmetry, 
the nonlinear atomic polarization should be parallel to 
the probe polarization (in the absence of magnetic field) . 
In consequence, all polarization components of the ND- 
FWM field are able to propagate along the PC beam. 

In order to allow the comparison between the theoreti- 
cal prediction and the experimental observation, we have 
incorporated into the calculation the intensity distribu- 
tion of the pump field in the atomic sample. Propagation 
effects were not considered. Such approach, expected to 
be valid for an optically thin sample, is only approxima- 
tive in our case since the peak probe absorption is around 
50%. Assuming that the F and B beams produce a per- 
fect standing wave, we have considered a sine wave distri- 
bution of Qi in the interval < f2i < ^imax- For each 
value of ill, the NDFWM field was calculated assuming 
that a steady state is reached. The total NDFWM field 
was taken as the sum of the contributions for each fli. Fi- 
nally the NDFWM power correspond to the square mod- 
ulus of the total field. The calculations were performed 
for an Fg = 4 ^ Fe = 5 transition with A = — 2r (corre- 
sponding to the experimental conditions) and 7 = O.OIF. 
No magnetic field was considered. The value of i^iMAX 
was adjusted to fit the experimental spectra. The best 
agreement was obtained for fliMAX = 18T. 

Figure |l2|b) represents the calculated spectrum for 
linear and parallel pump and probe polarizations with 
the parameters corresponding to the experimental condi- 
tions. The spectrum is dominated by the two sidebands. 



The width of the sidebands is mainly due to the inhomo- 
geneity of the Rabi frequency. The largest values of 
are responsible for the central peak presenting a width of 
the order of F while the atoms corresponding to fii < F 
are responsible for the narrow dip present at S — 0. The 
main features of the experimental spectrum are well re- 
produced. Some differences appear in the amplitude and 
shape of the narrow resonance around 6 = 0. According 
to the theoretical considerations presented above, this 
narrow resonance is the consequence of the coupling be- 
tween different Zeeman sublevels pairs through sponta- 
neous emission. Its width is governed by the time of 
flight relaxation rate 7 which effectively plays the role 
of a ground-state decay rate. For cold atoms, the aver- 
age time of flight across the exciting beams is rather long 
(more than I ms for a I mm diameter beam) and the cor- 
responding value of 7 < 3x I0~^ F too small to account 
for the observed resonance. However, for typical MOT 
temperatures, the average time os flight across one spa- 
tial period of the stationary wave produced by the F and 
B beams is three orders of magnitude shorter and corre- 
spond to 7 ~ 10~^ F as assumed in the calculation. At 
this point one should notice that a short travelling time 
across the standing wave pattern of the pump field in not 
compatible with our assumption of a steady state reached 
for each value of Oi. This suggests that a more sophisti- 
cated theoretical approach, incorporating the atomic mo- 
tion, would be more appropriate for a precise description 
of the spectra. In addition to the finite interaction time, 
another mechanism that can affect the narrow resonance 
around (5 = is the possible leakeage out of the closed 
two-level transition. In our case, this may occur though 
non resonant excitation of other excited state hyperfine 
levels. However the corresponding rate can be estimated 
to be smaller than IO~'^r. Also, calculations carried for 
open transitions give spectral shapes very different from 
those observed in Fig. |lj. 

The comparison between observed and calculated spec- 
tra in the case of linear and perpendicular pump and 
probe polarizations are presented in Fig. The main 
features of the spectra are well reproduced by the calcu- 
lation. Notice the significant increase in the maximum 
NDFWM yield with respect to the previous case. The 
spectrum is in the present case dominated by the narrow 
resonance at (5 = 0. The width of this resonance is also 
given at low pump intensities by 7. As already pointed in 
Q and [|l3|, this configuration provides the largest ND- 
FWM yield. The difference in the spectral profile existing 
between Figs. 
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and |I3| constitutes a clear demonstra- 
tion of the essential role of the Zeeman degeneracy and 
optical polarizations in the NDFWM process. 
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V. CONCLUSIONS. 

The process of NDFWM in a closed atomic transitions 
with Fe — Fg + 1 has been examined both theoreticaUy 
and experimcntaUy. Large differences in the spectra are 
observed for different choices of the exciting fields po- 
larizations revealing the crucial role of the internal level 
structure on the non-linear process. The spectra present 
distinct features which are determined by the different 
relaxation rates and characteristic frequencies of the sys- 
tem (F, 7, 1^1, A) H]. NDFWM was observed in a cold 
sample of cesium atoms in a PC experiment for two dif- 
ferent choices of the pump and probe polarizations. In 
spite of the difference existing between the experimen- 
tal conditions and the assumptions of the theory a good 
agreement between calculated and observed spectra was 
obtained. Nevertheless, some features of the spectra, as- 
sociated to the longest relaxation processes, indicate the 
need for a more detailed theoretical approach including 
the atomic motion and spatial field distribution. 
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FIG. 1. Schematic representation of the configuration con- 
sidered in the theory, e and g are degenerate levels. Solid 
vertical arrows represent the pump and probe fields (frequen- 
cies uji and iJ2 respectively). Dashed arrow: decay by sponta- 
neous emission. Oblique arrows: arrival and departure of the 
atoms into and from the interaction zone. The corresponding 
rates are indicated. 
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FIG. 2. Calculated NDFWM spectra for a pure two-level 
system, a) A = 0. b) A = 2r. 



FIG. 3. Different level configurations corresponding to the 
transition Fg = 1 — > = 2 for different choices of the pump 
and probe polarizations, a) Circular and equal, b) Linear 
parallel, c) Linear perpendicular. Solid (dashed) arrows cor- 
respond to the pump (probe) field. 
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FIG. 4. Calculated NDFWM spectra for the transition FIG. 5. Calculated NDFWM spectra for the transition 

Fg = 1 — > Fe = 2 for circular and equal pump and probe = 1 — > Fe = 2 for linear and parallel pump and probe 

polarizations (7 = O.OIF). a) A = 0. b) A = 2r. polarizations (7 = O.OIF). a) A = 0. b) A = 2F. 




FIG. 6. NDFWM power spectra for the transition 
Fg = I ^ Fe = 2 with A = 2r and 7 = O.OIF for differ- 
ent values of the pump field Rabi frequency Qi . 
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FIG. 8. Square modulus of the contribution to the atomic 
polarization at frequency ujs = uJi—5 from the pairs of Zeeman 
sublevels with \m\ = 1. 
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FIG. 10. Dressed atom energy level scheme for a transition 
_Fj = 1 — > _Fe = 2 driven by a tt polarized pump field. The 
notations are taken from Ref. The arrows illustrate two 
possible probe field transitions giving rise to resonances in the 
NDFWM spectra corresponding to opposite values of S. 
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FIG. 11. a) Experimental setup. AOM: acousto-optic mod- 
ulator. BS: beam splitter. PD: photodiode. MOT: mag- 
neto-optically trapped atomic sample. 



FIG. 12. Observed a) and calculated b) NDFWM spectra 
for linear and parallel pump and probe polarizations. The 
vertical axes scales are independent (the same scales are used 
inFigjl|). 
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FIG. 13. Observed a) and calculated b) NDFWM spectra 
for linear and perpendicular pump and probe polarizations. 
The vertical axes scales are independent (the same scales are 
used in Figjl2|). 
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